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Context. Radiative losses from optically thin plasma are an important ingredient for 
modeling plasma confined in the solar corona. Spectral models are continuously up- 
dated to include the emission from more spectral lines, with significant effects on ra- 



' diative losses, especially around 1 MK. 

Aims. We investigate the effect of changing the radiative losses temperature dependence 
due to upgrading of spectral codes on predictions obtained from modeling plasma con- 
■ fined in the solar corona. 

Methods. The hydrodynamic simulation of a pulse-heated loop strand is revisited com- 
paring results using an old and a recent radiative losses function. 

Results. We find significant changes in the plasma evolution during the late phases of 
plasma cooling: when the recent radiative loss curve is used, the plasma cooling rate in- 
creases significantly when temperatures reach 1-2 MK. Such more rapid cooling occurs 
when the plasma density is larger than a threshold value, and therefore in impulsive 
heating models that cause the loop plasma to become overdense. The fast cooling has 
the effect of steepening the slope of the emission measure distribution of coronal plas- 
mas with temperature at temperatures lower than ~ 2 MK. 

Conclusions. The effects of changes in the radiative losses curves can be important for 
modeling the late phases of the evolution of pulse-heated coronal loops, and, more in 
general, of thermally unstable optically thin plasmas. 

Key words. Methods: data analysis — Techniques: spectroscopic — Sun: corona — 
Sun: UV radiation — Sun: X-rays 



1 



Reale & Landi: Importance of radiative losses for coronal loops 
1. Introduction 



Coronal loops are the building blocks of the solar corona. They consist of pipe-like mag- 
netic field structures arching in the corona and connecting photospheric magnetic regions 
of opposite polarity. Loops are filled with optically thin, hot and relatively dense plasma, 
with temperatures ranging from 0.8 MK to a few million degrees, depending on the regions 
where they are located. There are two main classes of mechanisms that have been proposed 
to explai n the heating of th e loop plasma to coronal temperatures: steady (or "high fre- 
quency" , Warren et al.lboil ) heating mechanisms, and impulsive (or "low frequency" ) heat- 



ing mechanisms. Loop models that include steady heating predict a steady-state plasma, 
while in the pulse-heate d scenario the plasma evolves dynamically and spends most of the 



time in a cooling state (|Cargilllll994l ). In particular, a fast energy pulse might heat the 



plasma to more than 10 MK for a few seconds (e.g., 



Car gill 



2004; 



Reale k Orlando 



2008 



Guarrasi et al 



1994; 



Cargill k Klimchuk 



2010); then, the plasma is free to cool down 



losing energy by conduction to the chromosphere and by radiation. The balance between 
these two loss mechanisms is essentially determined by the local density: during and imme- 
diately after the heat pulse the loop plasma density is still low and conduction dominates, 
then denser plasma from the chromospher e fills the loop and radiation becomes the mos t 



efficient loss mechanism ( Antiochos 



1980; 



Cargill k Klimchuk 



2004; 



Reale 



2007 



2010) 



Observational e vidence seems to indicate that loop plasma might be heated impulsively 
(as reviewed in iKlimchukl 120061) . and recent studies using Hinode and Solar Dynamics 



Observatory data seem to co nfirm it ([Guarrasi et al 



2011 



Viall k Klimchuk 



2011). 



2010; 



Reale et al 



2011 



Terzo et al 



Regardless of the heating scenario, radiative losses are an important mechanism of 
plas ma cooling in coronal loops and need to be included with accuracy in loop models 



McClvmont k Canfield 



1983 



Antiochos k Sturrock 



J bood) 



1982; 



Bradshaw & Mason 



2003 



Mullcr et al. 2003; 



Bradshaw k Cargil 



are described as bundles of thin strands (Caigill 



a single short and intense heat pulse (jParker 



They become particularly important if loops 



1993 



Klimchuk 



2006), each ignited by 



19881 ) as the plasma cooling times are far 



longer than the expected duration of the heat pulse so that each strand spends most of 
its time cooling by radiation. Radiative losses are calculated by summing the emission of 
the plasma over the entire wavelength spectrum. The radiative emission of the optically 
thin plasma is dominated by the bremsstrahlung and free-bound recombination continua, 
and by line emission from all ions of all elements present in the coronal plasma. The total 
radiative losses are approximately proportional to the square of the electron density of th e 



plasma (|Tucker k 



Gould 



1966; 



Landini k Monsignori Fossi 



1970; 



Tucker k Koren 



1971) 



Landi k Landinil (|1999l ) showed that departures from such proportionality are lower than 
25%. The dependence on the electron temperature, on the contrary, is much more complex 
and, in typical coronal conditions, can be parameterized with a separate function of the 
temperature obtained by fitting a known curve to th e total radiative losses of isothermal 



spect ra calculated using a grid of temperature values ([Rosner et al 
199jJ). 



1978; 



Landi k Landini 
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The isothermal spectra used to determine the total radiative loss curve are obtained 
using a spectral code, which collects all the relevant atomic parameters and transition rates 
necessary to calculate the line and continuum emission of a n optically thin plasma. The 



1997 



Landi et al 



most popular co des available in the literature ar e CHIANT I ( Dcrc et al 
2012 ). AtomDB ( Foster et al1l2012 ). and SPEX ( Kaastra et al.lll996h : as new and improved 



calculations of atomic data and transition rates become available in the literature, these 
codes are costantly updated to extend their calculations to new lines neglected before, 
or improve the results for those already available. As a consequence, predicted radiative 
losses change, and differences from those calculated with earlier versions of the codes can 
sometimes be large. 

The consequences of upgrades to spectr al codes and radi ative losses can be impor- 
tant both for dat a analysis and diagnostics (|Testa et al.ll2012n . and for plasma modeling 



( Soler et al 



20121 ). For example, time-dependent hydrodynamic plasma models have been 



extensively used to describe loop plasma evolution. All models include the radiative losses 
in the energy balance, some times with a temper ature dependence des cribed with a piece 
wise p ower law function as in 



Rosner et al 



Raymond fc Smith 



1978), which is based on the 
( 1977t ) spectral model. More recent calculations of the radiative losses have a much larger 
loss rate at around 1 MK, due to higher metal abundances and the inclusion of a large 
amount of spectral lines from Fe vin-xv as well as other elements formed at similar tem- 
peratures. Differences at higher temperatures are more limited. Thus, for a plasma impul- 
sively heated to about 10 MK, we do not expect large differences for most of the evolution. 
However, when the plasma cools down and approaches 1 MK, the much larger loss rate 
will cause the loop to cool more efficiently. 

In this paper we will show in detail that, in a pulse-heated loop model, the enhanced 
radiative losses more easily and earlier lead to catastrophic cooling, that qualitatively 
changes the plasma evolution and may have important implications on both the observed 
emission distribution in active and quiet regions, and on the loop heating mechanisms 
as well. Our approach is to revisit a well-t ested multi-strand pulse-heated loop model 



(Reale & Orlando 



2008 



Guarrasi et al 



2010) using the standard and an updated radiative 



losses function and discuss the differences in the results. 

Section [5] summarizes the loop model focusing on the radiative losses function, Section[3] 
illustrates the results using different functions, Section [4] describes the implications for 
diagnostics, Section [5] discusses the results and their implications. 



2. The model 



Guarrasi et al 



(2010) that 



1982 



We consider a hydrodynamic simulation identical to the one in 
models a pulse-he ated loop strand with the Palermo-Harvard loop code (jPeres et al 
Betta et alll997 ). The loop strand is semicircular, symmetric with respect to the apex, and 
its symmetry axis is perpendicular to the solar surface. The loop half-length is L = 3 x 10 9 
cm and it includes a chromospheric layer at the footpoints that is linked to the corona 
through a steep transition region. 
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The plasma confined in each strand transports energy and moves only along the 



dependent hydrodynamic model (e.g., Nagai 


1980 




Peres et al. 


198z 


; Doschek et al. 


1982 




Nagai & Emslie 


1984; McClvmont & Canfield 




1983: 


MacNeicel 


1986; 


Gan et al. 


1991 


Hansteen 


1993 


: 


3etta et al. 1997; Antiochos et al. 


1999 


: Miiller et al. 2003: 


Bradshaw & Mason 2003 


Bradshaw & Cargill 


2006). through the equations ( 


Peres et al. 


1983: ] 


3etta et al. 


1997 


): 
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where n is the hydrogen number density; t is time, s is the field line coordinate; v is the 
plasma velocity; mu is the mass of hydrogen atom; p is the pressure; g is the component of 
gravity along the field line; fi is the effective coefficient of compressional viscosity (including 
numerical viscosity); f3 = n e /n is the ionization fraction where n e is the electron density; 
T is the temperature; k is the thermal conductivity (~ 9 x 10 -7 erg cm -1 s -1 K~ 7 / 2 ); ks 
is the Boltzmann constant; \ is the hydrogen ionization potential; P (T) is the radiative 
losses function per unit emission measure (discussed later); Q (s,t) is the power input per 
unit volume: 



'(M) =H + Hxf(t) 



(7) 



Ho is a low-regime (Ho = 3 x 10~ 5 erg cm~ 3 s^ 1 ) steady heating term which balances 
radiative and conductive losses for the static initial atmosphere. Hi is the amplitude of the 
heat pulse, and is assumed to be uniform along the loop. The time dependence of the heat 
pulse is a top-hat function, with f(t) = 1 for < t < 60s and f(t) = at any other time. 
The amplitude of the pulse is Hi = 0.38 erg cm -3 s _1 . We have checked that the presence 
of the steady heating term (Ho) is irrelevant for the entire strand evolution. For further 
investigation, we will also show simulations with different intensity of the heat pulse. 

The initial condition is that of a very low-pressure loop atmosphere, with a base pressure 
of po ~ 0.055 dyne cm -2 , which results in an apex temperature of To « 8.0 x 10 5 K. The 
chromosphere is assumed to be model F in 
is strictly maintained at all times. 



Vernazza et al 



(|1981l ). and its energy balance 



The Palermo-Harvard loop code (IPeres et al 



sively used for to mode l both flaring (IPeres et al.l 



loops (jReale et al 



2000; 



1982; 



Guarrasi et al 



3etta et al 



1987 



Betta et al 



1997) has been cxtcn- 



20011 ) and quiescent 



2010). The code has an adaptive mesh refinement 
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10 4 10 5 10 6 10 7 10 : 
T[K] 



Fig. 1. Radiative losses function P(T) (emissivity p er unit emission measure) as a func 



tion of temperature acco rding to 



Ros 



CHIANTI spectral code (jLandi et al 



neretal 



(|1978I ) (thin solid line) and to version 7 of 



20121 ) (thick solid line). 



(Bettaetal 



19971 ). to achieve adequately high resolution in the steep gradients along the 



strand and during the evolution. 

We replicate the simulation of 
functions P(T), sh own in Figure [U on e from 



Guarrasi et al 



(201 0j) wit h two d ifferent radiative losses 
|l978l) (hereafter RTV) com- 



Rosner et al 



puted according to IRavmon d & Smith (1977J), the other computed according to version 



120121) . assuming a density of 10 9 cm 3 and ioniza- 



Derd (I2009|). We have checked that radiative losses functions 
(|2008| ). fall in the range between these two 



Klimchuk et al 



7 of the CHIANTI code (|Landi et al 
tion equilibrium according to 
adopted in other works, e.g. 
curves and therefore we expect intermediate results when using them. 

The differences between the two curves are due to three sources. First, the atomic 
models in CHIAN T I are much larger and more sophisticated than those available to 



Raymond fe Smith! (jl977l) . and such larger atomic models provide vast amounts of ad- 



ditional spectral lines. The main differences fall around 0.5-3 MK, because of the increase 
in size of the Fe vni to Fexiv models, and around 10 MK, because of the larger models for 
Fe XVIII-XXIII. Increas es in the size of the mo dels for ions of other elements have a more 
limited effect. Overall. lLandi fe Landinil ([19991 ) found that increases in CHIANTI version 2 



atomic models and improvements on the atomic data caused radiative losses to change 
by as much as 50%. As even bigger models are used in CHIANTI version 7, we expect 
differences to be larger. 

The second source of variation lies in the ion fractions used in the two calculations. 
Ionization and recombination rates used to calculate the charge state composition at equi- 
librium are normally taken from theoretical calc ulations, which have been dramatically 
improved during the 30-40 years that separate the IRavmond fe Smith! (|l977l ) and the lat- 
est version of the CHIANTI code. The differences due to ion fr action improvements are 
expected to affect all elements and ions. lLandi fe Landinil (|l999f) . for example, estimated 
that different ion abundance datasets caused the radiative losses to change by up to 40%. 
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The third sou rce of variation are ele ment abundances. These are expected to pro- 
vide large effects: lLandi fe Landinil (|l999l ) reported differences of factors up to 2.5 when 
differ ent abundance da t asets are used. In the present comparison, the ab u ndanc es used 
by iRavmond fe Smith! (| 19T7T) are the cosmic values reported by lAllenl (|1973l) . while 
the CHIANTI code radiative losses were calculated using the coronal abundances by 
Feldman et al.l (]1992l ). The main difference between the two abundance datasets lies in the 
fact that the cosmic abundances of elements with First Ionization Potential (FIP) smaller 
than 10 eV have been increased in the coronal abundance dataset by a facto r »3.5 to 



account for the element fractionation in the solar corona known as the FIP effect (Feldman 



1992); additional, much smaller differences are found between the photospheric element 



abundances of ions with FIP > 10 eV and the cosmic values. The factor 3.5 enhancement 
involves most of the major elements emitting in the corona - Mg, Si and Fe - so that large 
effects are expected at temperatures in the 1-15 MK, where Fe emission dominates the 
spectrum. 



3. The results 



We have calculated the evolution of the loop-confined plasma for ~ 2000 s after 
the start of the impulsive heating. Such an evo l ution is known from many previous 



studies (e. 



Reale fe Orland 



Peres et al 



2008; 



1993 



Guarrasi et al 



Warren et al 
201ph 



2002 



2003 



Patsourakos fe Klimchuk 



2005 



Figure [2] shows samples of the temperature 
and density profiles along half of the strand at several different times (0s, 10s, 30s, 60s, 90s, 
120s, 300s, 800s and 1400s), that summarize the entire strand evolution. Figure [5] displays 
the simulation computed with the CHIANTI radiative losses, but the evolution obtained 
using the other loss curve is overall very similar. Because of the strong heating pulse, the 
plasma rapidly heats to above 10 MK (60s) along most of the loop. Then it gradually 
and uniformly cools down, and reaches ~ 0.1 MK, i.e. below the initial temperature, in 
about half an hour. The density follows a slightly different evolution and with a different 
timing. It increases initially with an evaporation front coming up from the chromosphere 
it = 10s, 30s). After the front has reached the loop apex, the density increases more uni- 
formly and reaches its maximum after about 5 minutes, i.e. much later than the end of the 
heat pulse. Then also the d ensity begins to decrease, due to draining driven by the cooling 



( Bradshaw fe Cargill 



20101 ) . After about half an hour, the density is still much higher than 



it was before the heat pulse. 

The effect of the different radiative loss functions is more significant at relatively late 
times. Such difference is illustrated in Figure[3J that shows the evolution of the temperature 
and density at the loop apex obtained using the two different loss functions. Both temper- 
ature and density at the loop apex are the same during the first ten minutes of evolution, 
when the temperature is larger than ~2 MK. After 10 minutes, the evolutionary tracks 
diverge: the plasma temperature and density calculated with the CHIANTI radiative losses 
decrease faster. 

The temperature evolution ends up with catastrophic cooling ( Parkerfl953 ; Field 1965 ) 



in both cases with temperature dropping abruptly by almost one order of magnitude to 
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0.5 1.0 1.5 2.0 2.5 
Coord, along the loop (10 9 cm) 



Fig. 2. Model temperature (top) and density (bottom) along half of the pulse-heated loop 
strand. The loop apex is at the right end of the X-axis. The plasma temperature and 
densities are displayed at several different times (in s) as reported in the two panels, from 
(heat pulse start) to 1400s. 



the minimum of ~ 2 x 10 4 K in little more than 3 minutes. Analytical approxi mations of 
the rate of decrease of th e plasma temperature have been derived in the past (jAntiochos 



1980; 



Cargill 



1994, 



1995!) and can be compared to those ob tained numer ically. Figure [3] 



shows two solutions obtained using the analytical expression (|Cargill 



1994): 



Ttop(t) — Tq 
where 



l/(l/2-«) 



3k B T 



(8) 



(9) 



n P(T) Q 

For both cases, we have assumed no = 7 x 10 9 cm -3 as density value , To = 6 MK as 
starting temperature of the radiative phase, from inspection of the simulations (see Fig. [3]), 
and a = —0.5 as effective power-law index of the radiative losses function approximate d 
with P(T) — \T a in the temperature range of interest, i.e. 5.5 < logT < 7 (Cargill 1994). 
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We have checked that the chosen values of no and Tq are in agreement with those obtained 



from equating the thermal conductive and radiative cooling timcscales (|Cargill fc Klimchuk 



2004). The two solutions that match well those of the simulations are obtained using two 
different values of the effective radiative losses, i.e. P(T)q — 1.5 x 10~ 22 and 1.1 x 10 -22 
erg cm 3 s -1 (see Fig. [T]). Substituting in Eq.© wc obtain r r fs 2400s and 3200s, for 
the faster (CHIANTI) and slower (RTV) cooling, respectively. The decay after the initial 
impulsive evolution, i.e. as soon as the temperature settles to about 6 MK, is generally well- 
described by Equation ©, although some details are lost. In particular, the catastrophic 
cooling that occurs at ~ 0.3 MK is well-reproduced. The numerical solution with CHIANTI 
predicts a faster decay below ~ 2 MK, that leads to catastrophic cooling at t ~ 1400 
s, i.e. ~ 600 s earlier than with RTV losses. This approximately corresponds to the two 
important changes of slope in the CHIANTI radiative loss function (Figured]) at T ~ 3 MK 
and 0.5 MK. Therefore, simulations are in good agreement with analytical descriptions, in 
which we change only the "effectiveness" of the radiative losses in a certain temperature 
range, i.e. average values of P(T). 

The density decreases more gradually than the temperature ( Bradshaw fc Cargilll2010 ). 



In the RTV solution we notice saw-toothed behaviour due to plasma bouncing back and 
forth along the closed loop during draining. For our purposes, it is more important to look 
at the average evolution. On average, the density obtained with CHIANTI losses decreases 
considerably faster than with RTV losses. The basic reason is that the more effective cooling 
makes pressure decrease faster in absolute value, and the pressure gradient, i.e. the force 
that sustains the plasma against draining by gravity, decreases faster as well, causing a 
faster draining. At this point the process is highly non-linear, i.e. catastrophic, and cannot 
be stopped even if the radiative cooling rate decreases with the density. Bradshaw and 
Cargill (2010) show the effect in terms of enthalpy, i.e. a peak of radiative losses is followed 
by a sharp increase of enthalpy losses. 

Overall, we might say that the density stays almost steady for more than 1500s 
with RTV losses, while its reduction becomes significant already after about 1000s with 
CHIANTI losses. This qualitative difference has important implications for the emission 
measure distributions. 

A very important qualitative difference of using CHIANTI instead of RTV radiative 
losses is the switching from regimes where there is a transition to catastrophic cooling 
to regimes where the rate of cooling remains roughly constant, i.e. both radiation and 
enthalpy losses increase smoothly with time, as shown in Bradshaw and Cargill (2010). 
This corresponds to the change of slope of the cooling occurring at about 2 MK. We 
have investigated in more detail what are the conditions required to trigger this change of 
behaviour. Figure 0] shows the evolution of the temperature and density at the loop apex 
obtained with the CHIANTI radiative losses and different values of the heat pulse intensity 
Hi. Higher Hi values lead both to higher apex maximum temperature and density values. 
Figure S] shows that the change of slope has a strong dependence on the density, regardless 
of the heat pulse intensity. In fact, the slope does not change if the loop plasma density 
remains approximately below the equilibrium density of a loop at a maximum temperature 
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500 1000 1500 2000 2500 

Time (s) 




1 000 1 500 

Time (s) 



2500 



Fig. 3. Evolution of the temperature (top) and density (bottom) at the lo op apex from 



Rosner et al 



hydro dynamic simulations using the two different radiative loss functions: 
(|1978r i (thin solid line) and version 7 of CHIANTI (thick solid line). Results obtained with 
analytical approximations (|Cargilllll994l ) are also shown (dotted lines, see text for details). 
The level of 2 MK is marked (dash-dotted horizontal line). 



of about 2 MK. According to iRosner et alJ (|1978f ) this density threshold can be obtained 
as: 



ng ~ — 
Lg 



(10) 



where ng is the density in units of 10 9 cm -3 and Lg is the loop half-length in units of 
10 9 cm. If the density is lower than this threshold value then the loop does not enter a 
catastrophic cooling phase. 
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Fig. 4. Evolution of the temperature (top) and density (bottom) at the loop apex obtained 
with the CHIANTI radiative losses and different values of the heat pulse intensity. The 
density equilibrium value (dashed horizontal line, bottom panel) for a loop at ~ 2 MK 
(dashed horizontal line, top panel) is marked. 



4. Diagnostic implications of different radiative loss curves 

The change of cooling rate has two important consequences. First, an important diagnostic 
result that can be obtained from coronal observations is the distribution of the emission 
measure as a function of temperature (hereafter EM(T)), 

rV(T) 

EM(T)= / n(TfdV (11) 
Jo 
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10000F 



1000 r 




1 I i I 

10 5 10 6 10 7 10 8 

T[K] 

Fig. 5. EM(T) curves obtained averaging the emission measure distribution vs temperature 
over the first 2000s of strand evolution, obtained from the simulation with the RTV (thin 
solid line) and the Chianti (thick solid line) radiative losses. The temperature of 2 MK is 
marked for reference (dotted vertical line). 



Typical coronal EM(T) distributions monotonically increase by a few orders of magnitude 
from ~ 0.1 to ~ 1 MK, have a broad peak around 2-3 MK and then decrease again to 
higher temperature , more steeply than they rise (e.g. 



Peres et al 



2000: 



Testa et al 



2011 



Warren et al 



201 lh . The EM(T) curve of a coronal loop can be predicted from theoretical 



models of loop strands like the one we used here by assuming that the loop consists of many 
strands, each heated randomly in time. Thus, the overall EM(T) curve will be determined 
by assuming that at any given time there will be loop strands at all stages of the strand 
evolution, and the total EM(T) curve can be obtained by averaging the EM(T) curve of a 
single strand over the entire strand evolution, in this case 2000 s. 

In Figure [S] we plot the space- and time-averaged EM(T) curves obtained from the 
simulations made using the RTV and the CHIANTI radiat i ve losses, in the cor onal part of 
the strand, i.e. excluding moss regions jGuarrasi et alJboioilwarren et alJhoill) . The figure 



shows that the high temperature portion of the EM(T) curve does not change in the two 
calculations, confirming that using different radiative loss curves does not change the results 
above 3-4 MK. On the contrary, the low temperature part is significantly different, due to 
the different slopes of the cooling at late times: the RTV simulation leads to a slightly higher 
peak and to a shallower distribution on the low temperature tail, that steepens below ~ 0.5 
MK. On the other hand, the faster cooling obtained with CHIANTI leads to a steeper curve 
below 2 MK, and thus to a sharper peak at 3 MK. The absolute difference between the two 
curves reaches a factor 5 around 1 MK; however, the difference in the slope of the two curves 
is even more important because such slope has been found to be an important parameter to 
discriminate betwe en low frequency and high frequency coronal heating mechanisms (e.g. 



Warren et al.ll2011l ). In the temperature range 5.9 < logT < 6.5 the slope in logarithmic 
1.1 for the RTV curve, and « 2.0 for the CHIANTI curve. These values are 



scale is 

smaller than those obtained by 



Warren et al 



(2011), but one cannot exclude that they 
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10000F 




10 5 10 6 10 7 10 8 

T[K] 

Fig. 6. EM(T) curves obtained averaging the emission measure distribution vs tempera- 
ture over the first 2000s of strand evolution, obtained from three simulations made using 
CHIANTI radiative losses and three different heat pulse amplitudes: the same used in 
Figure [5] (thick solid line), a lower one (lower thin solid line), and a larger one (upper thin 
solid line) (see also Fig. H]). 



may be compatible when all possible u ncertainties in data analysis and issues regarding 



DEM reconstruction (Testa et al 



20111 ) are taken into account. 



Moreover, the two different slopes on the cool side of the EM(T) curves lead to dif- 
ferent intensity rati os of EUV lines frequently observed by the available instrumentation 



(Wa rren et al 



201111 . Figure [5] shows two more distributions obtained with CHIANTI ra- 



diative losses and different heat pulse intensities, one with a smaller Hi value leading to 
a density steadily below the threshold value in Eq. (|10l) . one with a larger Hi leading to 
a very high density (Fig. H]). The figure clearly shows that the trend on the cool side of 
the EM(T) becomes shallower and more similar to the RTV results at low density: this 
means that i) the different slope is linked to the presence of different cooling rates, and 
ii) it can in principle be used to efficiently discriminate between coronal heating rates and 
mechanisms. 

The faster cooling with CHIANTI losses below 2 MK has a second, important impli- 
cation. The plasma in the strand spends a much longer time at temperatures well above 
1-2 MK than below them. If a loop is made by a multitude of strands, each heated once 
randomly, most of the strands within a loop will have a high temperatures at any given 
time, while only few of them can be observed at temperatures at or below ~ 1 MK. This 
might explain, at least in part, why active regions have been found to be mostly covered by 
hot ^3 MK loops with X-ray images and SDO/AIA 335 A channel images, and much less 
populated by warm -1 MK loops, like in the EUV 171 A channel in TRACE or SDO/AIA 
images, i.e. few loops are seen cooling from hot to warm status, much fewer than expected. 
To illustrate this effect, Figure [7] shows synthetic images of a multi-stranded loop in the 
two SDO channels obtained from the hydrodynamic simulations. We consider 200 straight 
strands grouped in bundles of 10 to mimic a limited instrument resolution. The straight as- 
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Fig. 7. Synthetic (straightened) loop images obtained from collecting groups of 10 ran- 
domly heated strands, in SDO/AIA 335A and 171 A channels, with models using RTV 
and CHIANTI radiative losses. The color scale is linear and exactly the same for a given 
channel (white is bright). The X-axis is measured in arcsec. 

pect may recall a loop projected on the solar disk. Each strand has the evolution described 
above, but the start time of the heat pulse is shuffled: the strands brighten at random 
times. We take a snapshot in a regime situation, i.e. the emission spatially averaged over 
the stranded loop is equal to the average emission over the whole strand time evolution. 
We show images obtained both with RTV and CHIANTI radiative losses. The color scale 
is linear and chosen to be exactly the same in each channel. In the 335A channel we see dif- 
ferences in the details but overall the loop appears to be bright. In the 17lA channel, there 
is a big qualitative change: the loop is much fainter, almost invisible, with the CHIANTI 
radiative losses. In a nutshell, the CHIANTI radiative losses significantly decrease the time 
that a given loop strand spends at a temperature between 1 and 3 MK, so that on the 
average it will be more difficult to observe loops in many narrow-band channels from AIA, 
TRACE, STEREO/EUVI and SOHO/EIT than expected when using RTV radiative losses. 
Since the onset of the accelerated (and then catastrophic) cooling depends on the plasma 
density, and hence on the rate of impulsive heating, changes in the radiative losses can 
have significant effects on heating rates determined from narrow-band images. 

5. Discussion and conclusions 

In this paper we have studied the effects of changes in the radiative loss curve on the 
evolution of an impulsively heated loop strand. Radiative losses can change due to 1) im- 
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provements in atomic models with the inclusion of more accurate atomic data and transition 
rates, or lines previously unavailable; 2) upgrades in the plasma ion abundance composi- 
tion; 3) changes in the plasma element abundance composition. We have carried out our 
tests using two radiative loss curves, whose differences were due to all three these causes. 
The most recent curve, from CHIANTI version 7, had much more efficient radiative losses 
at temperatures in the 0.5-3 MK range than the older one (RTV losses). We expect similar 
effects also using other losses curves obtained from other spectral codes, with similar atomic 
models and element abundances (SecU]). More efficient radiative losses may have effects on 
modeling other thermally unstable optically thin plasmas, such as flares, supernova rem- 
nants, accretion columns from circumstellar disks, novae, galactic cooling flows. Here we 
focus on their effects on the physics and structure of coronal loops, more in particular on 
the evolution of plasma confined inside a coronal magnetic flux tube and subject to a short 
and strong heat pulse. 

The effect is very small in the initial phases of the evolution. Most previous studies of 
the flaring plasma focused on the rise and initial decay phase of the flare, and much less on 
the very late phases, so their re sults are not affected by changes in the radiative losses. On 
the contrary. iReale et al.l (|2012| ) studied the EUV late phase of a flare using the older RTV 
radiative losses, and were able to reproduce the observed light curves with success. The 
agreement they found between model and observations suggests two possibilities. First, a 
small change of model parameters may be enough to maintain the same degree of agreement 
if the CHIANTI radiative loss function is used: for example, increased radiative losses 
may require a smaller heating to lead to the final catastrophic cooling that explained the 
observations. Second, the element abundances in the flaring plasma may be lower than 
the coronal values assumed in the CHIANTI rates used in the present work so that the 
differences between the RTV and the CHIANTI radiative losses will be smaller. 

Our results show that changes in the radiative loss curves have important effects on the 
late phases of a pulse-heated coronal loop strand, because they can change the cooling rate 
of the heated plasma and on the timing of the final catastrophic cooling. This can lead to 
immediate and interesting consequences on what we can expect from observations. 

First, the change of cooling rate occurs at relatively low temperature (~ 2 MK), and 
leads to considerable steepening of the cool side of the emission measure distribution. 
This change in EM(T) slope has important effects on our understanding of coronal plasma 
heating. In fact, earlier studies of pulse-heated loop models, using RTV-like radiative losses 
provided relatively flat EM(T) curves and such broad curves, not observed in active regions, 
have been invoked as an argument again st the effectiveness of impulsive heating in non- 
flaring active regions ( Warren et alj|2011 ). Thus, the change in slope caused by an improved 
radiative loss function can lead to a significant change of perspective. 

Second, the temperature of the cooling plasma decreases more rapidly below 2 MK, 
and this might explain why we see very few loops cooling from X-rays to EUV bands in 
active regions. 

Third, it is also interesting to notice that we do not expect faster cooling to occur in 
lower-regime loops: if the loop is heated by less intense pulses or in a more gradual and 
longer-lasting way, the plasma might never start a fast cooling. Whenever we observe 1 
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MK loops in UV "warm" channels this might provide an upper limit to the intensity and 
duration of the heat pu l ses th at provide the energy. 



Fourth, 



Reale et al 



(|2012f) showed that a catastrophic cooling occurred at the end of 
the evolution of a post-flare loo p, and explained the so-called EUV late phase observed by 
EVE during the decay of flares (jWoods et al.ll201ll) . Changes in radiative losses cause the 
predictions of the EUV late phase of flares to change considerably, and affect any estimate 
of flare impulsive heating rates made combining cooling loop models with observed EUV 
light curves. 

This work uses one the latest release of an up-to-date spectral model. However, chances 
are that spectral codes will keep evolving and include many more lines and more accurate 
data as they become available . For example, th ere is evidence that we may still miss 



many important spectral lines ([Testa et al 



20121 ) that can have significant effects on the 



radiative losses. As a consequence, we expect even more significant effects due to further 
spectral model improvements, so the present work can even be considered conservative to 
this respect. Also, element abundances are very important for radiative losses, so studies 
that provide better constraints on abundances are highly encouraged. On the contrary, 
effects of deviations from ionization equilibrium of the emitting species are negligible in 
the late phase of the flare, because the plasma is still at relatively high density, so that the 
ionization equilibrium times are quite short and the plasma is abl e to respond quickly and 



adapt to the rapid changes in the temperature (jReale et al 



2012). 



In conclusion, in this work we report on the effect of upgrading the spectral models of 
optically thin emitting plasmas on modeling coronal plasma. In particular, we find that, 
while the increased emissivity obtained with new models has negligible effects on most of 
the plasma evolution predicted by models and in many plasma conditions, there is a very 
important implication for modeling of plasma confined in coronal loop strands heated by 
strong and fast energy pulses. Enhanced plasma emission in the 0.5-3 MK temperature 
range leads to faster cooling that might explain several pieces of evidence that questioned 
the validity of pulse-heated loop models. We point out that this kind of effects might 
be important also for modeling other thermally unstable optically thin plasmas, such as 
supernova remnants, accretion columns from circumstellar disks, novae. 
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